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Abstract. The spatial sensor characteristics of a 6cm TDR flat band cable sensor 
section was simulated with finite element modelling (High Frequency Structure 
Simulator-HFSS) under certain conditions: (i) in direct contact to the surrounding 
material (air, water of different salinities, different synthetic and natural soils (sand- 
silt-clay mixtures)), (ii) with consideration of a defined gap of different size filled 
with air or water and (iii) the cable sensor pressed at a borehole-wall. The 
complex dielectric permittivity e*{uj,Ti) or complex electrical conductivity a*(uj,Ti) = 
iu!e*{uj,Ti) of the investigated saturated and unsaturated soils was examined in the 
frequency range 50MHz-20GHz at room temperature and atmospheric pressure with 
a HP8720D- network analyser. Three soil-specific relaxation processes are assumed to 
act in the investigated frequency-temperature-pressure range: one primary a-process 
(main water relaxation) and two secondary (a', /3)-processes due to clay- water- ion 
interactions (bound water relaxation and the Maxwell- Wagner effect). The dielectric 
relaxation behaviour of every process is described with the use of a simple fractional 
relaxation model. 3D finite element simulation is performed with a A/3 based adaptive 
mesh refinement at a solution frequency of IMHz, lOMHz, O.lGHz, IGHz and 12.5GHz. 
The electromagnetic field distribution, S-parameter and step responses were examined. 
The simulation adequately reproduces the spatial and temporal electrical and magnetic 
field distribution. High-lossy soils cause depending on increasing gravimetric water 
content and bulk density an increase of TDR signal rise time as well as a strong 
absorption of multiple reflections. Air or water gap work as quasi wave-guide, i.e. the 
influence by surrounding medium is strongly reduced. Appropriate TDR-travel-time 
distortions can be quantified. 

Keywords: lossy dielectrics, finite element modelling, HFSS, dielectric spectroscopy, 
fractional relaxation 
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1. Introduction 

Soil science, geophysical prospecting, agriculture, hydrology, archeology and 
geotechnical engineering have benefited greatly from developments in radio and 
microwave technology. Electromagnetic techniques are used to estimate soil and rock 
physical characteristics such as water content, density and porosity ([17], [37], [27], 
[25]). Both invasive methods, such as time domain reflectometry ([19], [IQ], [15]) and 
cross borehole radar [in], and noninvasive methods, such as capacity methods ([IS], 
[25] , [27] ) and ground penetrating radar ([21], [2], [12], [6], [13]) are used. Common 
to all these techniques is the fact that electromagnetic wave interaction depends on 
dielectric properties of rock or soil deposit through which it travels, which are influenced 
by chemical composition, mineralogy, structure, porosity, geological age and forming 
conditions. Besides, several additions like ubiquitous water have an effect on the 
dielectric properties. 

In particular, knowledge of the spatial and temporal variability of water saturation 
in soils is important to obtain improved estimates of water flow (and its dissolved 
components) through the vadose zone. Due to its accuracy and potential for automated 
measurement, TDR has become one of the standard methods to measure spatial and 
temporal variability of water contents in laboratory soil cores and experimental field 
plots [15j|. For this purpose, the object of numerous experimental and theoretical 
investigations is the development of general dielectric mixing models for a broad class 
of soil textures and structures ([IS], [13], ^0]). Mostly, these empirical, numerical or 
theoretical models base on the assumption of a constant dielectric permittivity of the 
soil as a function of volumetric water content in a narrow frequency range around IGHz 
([S], [IH], [1], [9], [39]). However, the strong frequency dependence in the dielectric 
relaxation behaviour below IGHz due to a certain amount of swelling clay minerals in 
nearly each real soil is considered only insufficiently ([30], [15], [32], [25]). 

The type of multi-scale structure renders the analysis of dynamic data in clays 
rather complex. The problem has been addressed both by experimental and modelling 
techniques. Besides broadband dielectric spectroscopy on clay-water suspensions 
microscopic simulations of clays have been an active field of research since the late 1980s 
and began with simulations at ambient temperature and pressure, of clays with various 
cationic species. More recently, several studies appeared dealing with non-ambient 
conditions (increased temperatures and pressures), which are primarily linked to the 
issue of storage of radioactive waste or bore- hole stability [35]. Previous experimental 
and modelling results suggest that clay-water systems have multiple relaxation processes, 
such as interfacial polarizations around the clay particles and rotational relaxation of 
bound and free H2O. Therefore, the dielectric behavior is expected to be complicated. 
Useful and precise dielectric information may only be obtained when each relaxation 
process is extracted from the complicated overall behavior based on the measurement of 
the complex dielectric permittivity over a broad frequency range and at high resolutions 
([22], [21], [31], [33]). 



Determination of the spatial TDR-sensor characteristics. 



3 



In the present study, the the dielectric relaxation behaviour of the investigated 
saturated and unsaturated soils was examined in the frequency range 50MHz-20GHz. To 
parameterize the dielectric spectra three soil-specific relaxation processes are assumed 
to act as a function of water saturation and porosity in the investigated frequency- 
temperature-pressure range: one primary a-process (main water relaxation) and two 
secondary (a', /5)-processes due to clay-water-ion interactions (bound water relaxation 
and the Maxwell- Wagner effect). The dielectric relaxation behaviour of every process is 
described with the use of a simple fractional relaxation model ([22], [21], [23], [I9], [T6]). 
The chosen approach enables a characterization of the dielectric relaxation behaviour 
with the separation in the observed relaxation processes in dependence of the porosity 
and the water content. 

A simplification frequently utilized in TDR applications is the use of an idealized 
equivalent circuit for the sensor without consideration of losses due to the skin-effect 
or radiation from the sensor as well as the assumption of a homogeneous sensitivity 
distribution along the sensor ([15], [20], [15], [2H])- In addition, a frequently arising 
problem in various applications is the direct contact between sensor, e.g. a flexible 
band cable, and surrounding medium. An air or water gap between sensor and soil 
leads to dramatic under or overestimation of water content. A suitable tool for an 
examination of this specific problems offers three-dimensional numeric finite element 
simulation ([28]). For these reasons, in this study the spatial sensor characteristics of 
a 6cm TDR flat band cable sensor section was simulated with electromagnetic finite 
element modelling (Ansoft-HFSS"^*"'^ , High Frequency Structure Simulator). In order 
to carry out the finite element calculations as realistically as possible, the measured 
frequency-dependent dielectric permittivity was considered. Moreover, the simulations 
were performed under certain conditions: (i) in direct contact to surrounding material, 
(ii) with consideration of a defined gap of variable size filled with air or water and (iii) 
cable sensor pressed at a borehole-wall. 



2. Theoretical Background 

Time domain refiectometry measures the propagation velocity of a broadband step 
voltage pulse (typical values: rise time t^. ^ 70ps, sampling increment At ^ 20ps) with 
a bandwidth of around 20kHz to 25GHz (Nyquist-frequency: fmax = 0.5/At). But due 
to the limitations of used connectors, type of the TDR device, coaxial cable type and 
length the effective bandwidth is reduced distinctly [31]. Under atmospheric conditions 
the velocity of this signal is a function of the frequency uj = 27r/ and temperature T 
dependent effective relative complex permittivity ^^f^i^, T) = ^'^f^i^, T) ~ J^gff('^, T) of 
the material through which it travels. The overall losses £:^g(u;,T) = e'^{uj,T) + '^^^^^^ 
of the material which have to be considered result from the dielectric losses £^'(co',T) 
and the conductive losses due to a direct current electrical conductivity aj:,c{T). 

Here, Eq is the permittivity of free space and = —1 is the imaginary unit. It is often 
convenient to consider the analogy of propagation phase velocity and attenuation of an 
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electromagnetic plane wave (see [IH], [S], [U]) IS])- 

vp{u, T) = cV2 ^Je'^^{uJ, Tf + e^^iu, Tf + e'^^{uj, T)^ (1) 

/?(^, T) = uj^e'^^iu, Tf + e;;g(a;TT)^ - T) ■ [cV^f' (2) 

where c = ^Jsq^q'^ is the velocity of light and /xq the magnetic permeability of vacuum. 
Hence, any modulation of an electromagnetic wave in a real medium will 
a group velocity according to the Rayleigh equation ([H], [12]): 

(3) 



dk 



Vg = — = Vp 



f dvp 



Vp df 

herein, k denotes the wave number. The flat band cable of length / consists of three strip 
conductors embedded in a polyethylene band (see [20], [27] for details). The effective 
group or phase velocity of the signal v'p in a perfect dielectric (pure real dielectric 
constant Er = e*g = const without dispersion and conducting losses) surrounding the 
cable sensor is in principle only a crude approximation of the real material properties 
especially at frequencies / < IGHz and / > lOGHz (c.f. Fig. [1]) 

21 c 



v'p = V 



9 



t 



(4) 



where t is two way travel time. Considering anomalous dispersion equation (jl]) is referred 
to as a high frequency approximation of phase velocity (Fig. [1]). In contrast the high 
frequency attenuation approximation j3h {oj^T) for real soils frequently used in ground 
penetrating radar applications works considerably well (c.f. [6j, Fig. [1]) 

with impedance Zq = cfiQ of vacuum. We now consider the soil as a four-phase medium 
composed of: air, quartz grain, water and clay. In the particular case of spatial TDR 
the surrounding medium in the direction oft the band cable is described by a relative 
effective permittivity £*g(x, a;,T, p). It depends on position x, angular frequency u and 
contribution due to several relaxation processes via relaxation time t{T,p) on absolute 
temperature T and pressure p ([22], [30]) 



Ti{T,p) = Ki——exp — - — (6) 



kpT ^ \ RT 

Herein, h denotes the Planck- const ant, fc^-Boltzmann constant, ^ 1 the transmission 
coeffitient, R gas constant and Ea^i{T,p) = AGi{T,p) + TASi{T,p) activation energy 
with free enthalpy AGj(T,p) and activation entropy ASi{T,p) of the z-th process (j22j. 
[2T]). Dielectric loss spectra of saturated and unsaturated soils very often show a marked 
deviation from simple Debye-behaviour ([H], [Hj, [22], [25j). Based on the theory 
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Figure 1. (top) Phase velocity wp, liigh frequency approximation u^, corresponding 
group velocity Vg as well as (bottom) attenuation (3 and corresponding high frequency 
approximation [3h of (left) a sand-bentonite mixture (SB 50/50 1-6) with different 
gravimetric water content 8g and bulk density g (see section [31 Tab. [31 Fig. [S|). (right) 
Comparison of natural water with aoc — 0.18S'/m, a mikrosil quartz powder and SB 
50/50-4 (see Tab. [Hand [3] below). The calculations were carried out under utilization 
of the measured dielectric permittivity. 



of fractional time evolutions Hilfer [16] derived a Jonscher type function ^23j for the 
complex frequency dependent dielectric permittivity of amorphous and glassy materials 

^eff,«(^' - = T- — ^^7T7^ — ^ 

with high frequency limit of permittivity £00 ? relaxation strength Aei{T) as a function 
of temperature, angular frequency oj and stretching exponents < ai, [3i similar to 
the familiar empirical Havriliak and Negami [l4j, Cole-Cole [3], Cole-Davidson [5] or 
Kohlrausch- Williams- Watts ([26], [51]) dispersion and absorption functions. For the 
particular case = and A = 1 (I7j) transforms to the Debye model. 



3. Material and Experiments 

The complex effective dielectric permittivity e*g of saturated and unsaturated soils was 
examined in the frequency range 50MHz-20GHz at room temperature and atmospheric 
pressure with a HP8720D- network analyser. This was performed using a combination 
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Figure 2. Schematic diagram of the experimental set-up for measuring the dielectric 
permittivity: (top) longitudinal section of the coaxial transmission line with an outer 
diameter of the inner conductor di = 7mm, an inner diameter of the outer conductor 
do — 16mm and a length / — lOOmm, (bottom) open-ended coaxial-line (HP85070B). 



of open-ended coaxial-line (HP85070B) and coaxial transmission line technique (sample 
holder (7xl6xl00)mm^) (Fig. [2]). The synthetic soil samples were incrementally wetted 
from air dry up to saturation with natural water and equilibrated 12h. From the 
prepared sample a subsample was taken with a retaining ring (hight 100 mm, inner 
diameter 100 mm). The retaining ring was used as the sample holder for the HP85070B 
probe. Care was taken to pack the soil in the transmission line to a homogeneous bulk 
density g and to a constant volume. After each dielectric measurement bulk density 
g as well as gravimetric water content dg were determined. The measured complex 
S-parameter values Sij were used to calculate complex dielectric permittivity with a 
commercial software (HP 85070/71C Materials Measurement Software) after calibration 
with open, short, and load standards. 

Different natural and synthetic soils were investigated. Here, we present our results 
for synthetic soil SB50/50. It is a mixture of 50wt.% sand (grain size <2mm) and 50 
wt.% bentonite (Calcigel: 71wt.% Ca- dioctahedral smectite, 9wt.% illite/dioctahedral 
mica, lwt.%kaoline, lwt.% chlorite, 9wt.% quartz, 5wt.% feldspar, 2wt.% calcite, 2wt.% 
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Figure 3. SEM-EDX micrograph of tlie investigated bentonite Calcigel, Philips XL30 
ESEM-FEG. 

Table 1. Physical properties of the investigated bentonite 8 ) in comparison to the 
commercial Mikrosil quartz powder (type 350, median grain size 11/im). 





Ca-Bentonite 


Mikrosil 


Dry density Qd 
Specific surface area 
Density 

Cation excliange capacity 


1.33 g/cw? 
493 mVg('') 
2.847 g/cm^ 
62.0 meq/100g(^) 


0.69 g/cm^ 
0.38 mVg 
2.65 g/cm^ 



dolomite). 

Ca- dioctaliedral smectites are clay minerals, i.e. they consist of individual 
crystallites the majority of which are <2/im in largest dimension. The crystal structure 
of clays has been established from X-ray diffraction studies for almost all types of 
common clays. Smectite crystallites themselves are three-layer clay minerals. Individual 
clay layers consist of fused sheets of octahedra of AF"*" or Mg^"*" oxides and tetrahedra 
of Si^+ oxides. Substitution of k\^+ , Mg^+ or Si^ with lower valence ions results in 
an overall negative charge on the layers, which is then compensated by cationic species 
(counterions, in the current case Ca^^) between clay layers (interlayers) ([29], Fig. [3]). 

Under increasing relative humidity the compensating ions become hydrated and the 
spacing between individual layers increases. While the detailed swelling characteristics 
of a smectite (interlayer spacing as a function of relative humidity) depends crucially 
on the charge of the clay layers (magnitude and localisation) and the nature of the 
compensating ion, in general, it occurs in three stages. Swelling begins in a step- 
wise manner (discrete layers of water formed in the interlayer, states referred to as 
monolayer/monohydrated, bilayer/ bihydrated, etc.), becomes continuous thereafter and 
in the extreme a colloidal suspension of clay particles (<10 aligned layers) is formed. 
Beyond the microscopic scale, aggregates of aligned clay layers form particles of the 
order of 10 nm - 1000 nm in size, with porosities on the meso- (8 nm - 60 nm) and 
macroscale (>60 nm) ([35j and citations in it). In Tab. [1] the physical and chemical 
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Figure 4. Complex relative dielectric permittivity e*g: = — -^^eff ^ function 
of frequency of sand-bentonite-mixture SB50/50-4 with GDR-fit under consideration 
of three relaxation processes: one primary a-process (main water relaxation) and two 
secondary (a', /?)-processes due to clay- water- ion interactions (bound water relaxation 
and the Maxwell- Wagner effect). In this way it was possible to separate the singe 
relaxation processes and to extrapolate the dielectric permittivity to IMHz. 



properties of the used bentonite calcigel are summarized. 

Three relaxation processes are assumed to act in the investigated frequency- 
temperature-pressure range: one primary a-process (main water relaxation) and two 
secondary (a', /3) -pro cesses due to clay-water-ion interactions (bound water relaxation 
and the Maxwell- Wagner effect, see Fig. Hj). The effective permittivity of a multiphase 
soil-mixture can be determined by the complex relative permittivity of water , 
bound water e^^, the contribution due to clay- water- ion interaction e*^^^ as well as the 
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Figure 5. Dispersion and absorption curve of sand-bentonite-mixture (SB50/50) for 
seven gravimetric water contents Og and bulk densities q. 



real and constant permittivity of quartz grain Esand and air ([47], [37], [44], [48], [22] . 
|21] . [1]). The dielectric relaxation behaviour of each process is described by a fractional 
relaxation model according to ([7]) considering relaxation time distributions H{t). This 
allows the complete spectrum to fit as a function of water content 6g and bulk density 
g at constant temperature and pressure with the use of a generalized dielectric response 
(GDR, Fig. a and [5]): 

3 



A^i , a DC 

A Shuffled Complex Evolution Metropolis (SCEM-UA) algorithm is used to find 



e:d^)-e^ = y ^ J- -3 -^ (8) 
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Table 2. Parameters of the three relaxation processes from GDR- fitting (i = 
[a, a',/?]); gravimetric water content w, bulk density p, porosity n, water saturation 
Sw and volumetric water content 9 = S^n of SB-50/50. 





50-0 


50-1 


50-2 


50-3 


50-4 


50-5 


50-6 


Mikrosil 


w [%] 


0.6 


5 


11.4 


18.1 


28.1 


33.5 


41.5 


25.47 




1.35 


1.14 


1.37 


1.34 


1.75 


1.52 


1.66 


1.55 


71 


0.5 


0.6 


0.55 


0.59 


0.53 


0.62 


0.64 


0.56 




0.02 


0.09 


0.25 


0.34 


0.66 


0.54 


0.63 


0.71 


9 \%] 


0.81 


5.42 


13.84 


19.86 


35.36 


33.86 


40.3 


39.48 


^ oo 


1.03 


1.15 


1.23 


2.07 


3.38 


3.77 


6.69 


1.26 


Af 


0.86 


0.05 


0.01 


0.13 


5.44 


3.67 


19.57 


16.73 


Ivs] 

' CK J 


4.99 


9.18 


8.44 


7.43 


6.06 


8.74 


9.8 


9.25 




0.004 


0.004 


0.07 


0.025 


0.075 


0.091 


0.095 


0.002 


aa (fixed) 




























0.25 


2.16 


3.18 


5.06 


10.31 


19.7 


5.77 


2.09 




5.37 


9.64 


9.34 


9.63 


9.93 


9.7 


9.32 


0.27 




0.004 


0.081 


0.081 


0.075 


0.003 


0.082 


0.034 


0.048 


Q!q,' (fixed) 




























0.05 


1.57 


4.9 


11.06 


43.35 


62.17 


74.99 


0.87 


T/3 [ns] 


98.46 


0.5 


0.6 


0.51 


0.52 


0.65 


0.83 


22.6 


1-/3/3 


0.65 


0.99 


0.56 


0.69 


0.66 


0.68 


0.74 


0.92 


a/3 (fixed) 


1 


1 


1 


1 


1 


1 


1 


1 


CDC [mS/cm] 


6.68E-5 


0.03 


0.05 


0.75 


3.37 


4.62 


5.31 


0.12 



best GDR fitting parameters [I5], Tab. [3]). This algorithm is an adaptive evolutionary 
Monte Carlo Markov Chain method inspired by the SCE-UA global optimization 
algorithm of [7j and combines the strengths of the Metropolis algorithm [36j , controlled 
random search [38], competitive evolution [18], and complex shuffling p'J to obtain 
an efficient estimate of the most optimal parameter set, and its underlying posterior 
distribution, within a single optimization run. The resulting relative error of each 
parameter is less than 3%. 

The relaxation parameters obtained from GDR-fit are presented in Fig. El 
Relaxation strength Asi of each process and relative high frequency permittivity Soo 
depend on moisture content. Above a gravimetric water content of ~ 30wt.% the a'- 
process decrease and the a-process strongly increase. This suggest an increase of the 
primary a-relaxation at the expense of the bound water process a'. In contrast at the 
highest saturation level of 0.66 at a porosity of 0.53 primary a and low frequency (3- 
process decrease. Relaxation time Tj and distribution parameter /3i are nearly constant. 
An exception represents sample 50-0 which shows clear deviation from the general trend 
but because of the very low gravimetric water content of 0.6 wt%. Without swelling 
an increasing saturation reduces the gas-filled pore space, while the effective pore space 
available for the fluid becomes larger. This leads to the simple calibration models 
mentioned in the introduction to determine the volumetric water content with dielectric 
measurements. Due to the swelling of the clay minerals in the process of the hydration 
the distribution of immobile and effective pore space cannot be considered as constant. A 
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Figure 6. (top) Relaxation strength Ae^ of the i-th process in comparison to the 
relative high frequency permittivity function (left) of gravimetric water 

content and (right) water saturation Sw (bottom) Apparent direct current electrical 
conductivity <jdc- 



circumstance which comphcate a careful reahzation of accurate measurements is the fact 
that a multitude of variables of the measurement conditions and the sample preparation 
essential determine the intensity of the swelling and the resulting change of the pore 
structure ([Ij). 

The results show the potential of the chosen approach but a detailed explanation 
of this complex behavior is beyond the scope of this paper. In general, there is a need of 
further systematic investigations by broadband dielectric spectroscopy of saturated and 
unsaturated soils under consideration of the swelling process and with an utilisation of 
microscopic modelling. 



4. HFSS-Simulations 

The transfer or scattering function Sij{uj) of the flat band cable section (Fig. [7]) 
was simulated by finite element modelling (commercial software from Ansoft: High 
Frequency Structure Simulator-HFS^ under certain conditions: (i) in direct contact to 
the surrounding material (air, water of various salinities, different synthetic and natural 
soils (sand- silt-clay mixtures)), (ii) with consideration of a defined gap of various size 
(total high 2 mm, 3 mm, 5 mm, 7 mm or 10 mm) filled with air or distilled water and 
(iii) cable sensor pressed at a borehole-wall. 

In HFSS tangential element basis function interpolates field values from both nodal 
values at vertices and on edges. Surfaces of the structure (air-box) in the y-z an x-y 



X HFSS is a standard simulation package for electromagnetic design and optimization. 
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plane are radiation boundaries and the second-order radiation boundary condition is 
used 

V X = jkoEt - f Vi X (Vt X Et) + f Vt(Vt • Et) (9) 

where Et is the component of the E-ficld that is tangential to the surface and /cq is 
the free space phase constant. The second-order radiation boundary condition is an 
approximation of free space. The accuracy of the approximation depends on the distance 
between the boundary and the object from which the radiation emanates. For this reason 
a sensitivity analysis was carried out for a total height of the structure between 20 mm 
and 50 mm in dependence of the material properties. The influence of the boundary 





Figure 7. Model geometry of (top) flat band cable surrounded by saturated and 
unsaturated soil with a gap filled with air or water and (bottom) cable sensor pressed 
at a borehole-wall with air or soil as borehole filling. 
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layer for the simulation into air can be neglected for a distance from the cable sensor 
greater than 12 mm so a minimum height of 25 mm is used. 

Surfaces in the x-z plane are wave ports. HFSS assumes that each wave port 
is connected to a semi-infinitely long waveguide that has the same cross-section and 
material properties as the port. When solving for the S-parameters, HFSS assumes 
that the structure is excited by the natural field patterns (modes) associated with these 
cross-sections. The 2D field solutions generated for each wave port serve as boundary 
conditions at those ports for the 3D problem. The final field solution computed must 
match the 2D field pattern at each port. 




Figure 8. Electric field distribution @\2.bGH z for the investigated flat band cable 
surrounded by air, sand-bentonite-mixture (SB 50/50-4) with 6g =28.14wt.% and 
p = l.Tdg/cm^ as well as a defined 3mm air or water gap. (left) Cross section, (right) 
longitudinal section of middle conductor. 



The simulation is performed with a A/3 based adaptive mesh refinement at solution 
frequency of 1 MHz, 10 MHz, 0.1 GHz, 1 GHz and 12.5 GHz. Broadband complex S- 
Parameter are calculated with an interpolating sweep in frequency range 1 MHz - 12.5 
GHz with extrapolation to DC. The electromagnetic field distribution, S-parameter 
and step response (200 ps rise time) of the structure were computed in reflection and 
transmission mode. 



5. Discussion 



The simulation adequately reproduces the spatial and temporal electrical and magnetic 
field distribution in comparison with 2D-FE investigations of [20] (Fig. [8l). Fig. [9] 
represents reflection and transmission factor as well as the corresponding TDR waveform 
for the sand-bentonite mixture at various water contents and bulk densities in reflection 
and transmission mode. As a reference material the simulation results for the cable 
sensor surrounded by air are included. As expected, the appropriate resonances in the 
reflection coefficient Su shift with rising water content to deeper frequencies and the 
attenuation increases. In time domain onset travel-time as well as rise time of the TDR 
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Figure 9. (left) Input return loss magnitude or reflection coefficient 5*11 and forward 
transmission or transmission coefficient as well as (right) TDR-waveform in 
reflection and transmission mode for simulated flat band cable structure, surrounded 
by air and sand bentonite mixture of various water contents and bulk densities (see 
Tab. El). 



signals are analyzed. As is clearly recognizable in both reflection and transmission mode, 
onset time increase and rise time decrease with rising water content. 

Qualitatively the numerical calculation shows that sensitivity characteristic of the 
cable sensor changes along the sensor in dependents of the dielectric relaxation behavior 
of the surrounding material (Fig. [8]). The investigated high- lossy sand-bentonite 
mixture cause in dependence of increasing gravimetric water content 9g and bulk density 
Q an increase of TDR signal rise time as well as a strong absorption of multiple reflections. 
This leads to a frequency dependent decrease of spatial resolution and penetration depth 
(sensitivity region around the sensor) along the fiat band cable and fixed maximal length 
of the moisture-sensor available for application. 

Coupling problems caused by air or water gaps lead to dramatic travel time 
distortion even for very small gaps (Fig. [TOl) . The air filled gap with a thickness of 
0.25mm on both sides of the cable sensor already leads to the drastic underestimation 
of water content of 36%. In contrast a drastic overestimation occurs in the case of a 
water filled gap for the same gap size. Further increasing gap size leads to a maximum 
in the characteristic (anomal behavior) at a gap thicknesses of approximately 2.7 mm, 
whereby up to the complete filling of the structure with water the effective dielectric 
constant decrease and the rise time increase. An explanation of this effect and the 
consequences linked with it are however still pending and require further numerical and 
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experimental investigations. Moreover, with HFSS it is not possible to consider an 
exchange of material between air or water-filled gap and surrounding medium. 

As a consequence the gap work as quasi waveguide, i.e. the infiuence of the 
surrounding medium is strongly reduced but changes in dielectric properties along the 
cable sensor are reproduced. 
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Figure 10. Real effective relative permittivity Egg- = SZ/^Q^ggi; according to equation 
(|4]) plotted against TDR rise time i^jg^ (in reflection mode) for all sensor configurations 
and investigated cases. 



6. Conclusion 

In this study, the spatial sensor characteristics of a 6cm TDR fiat band cable section is 
simulated with finite element modelling in combination with dielectric spectroscopy. 

For this purpose the dielectric relaxation behaviour of saturated and unsaturated 
soils is examined in the frequency range 50MHz-20GHz. The dielectric relaxation 
behaviour is described with the use of a generalized fractional relaxation model. With 
this approach the frequency depended dielectric permittivity is determined based on a 
parametrization of each relaxation processes as a function of water content and porosity. 
This enables a development of improved calibration strategies. However, there is a need 
of further systematic investigations by broadband dielectric spectroscopy of saturated 
and unsaturated soils under consideration of the swelling process and with an utilisation 
of microscopic modelling. 

The three-dimensional numeric finite element simulation in HFSS provide an 
informative basis of the sensor characteristics under consideration of the frequency 
dependence of the measured complex dielectric permittivity. It is shown, that an air 
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or water gap between sensor and soil leads to dramatic under or overestimation of 
water content already for a gap thickness of 0.25 mm upper and below the cable sensor. 
Therefore, the application of the flat cable as a moisture sensor requires an accurate 
installation technique of cable-like elements (especially for long sensors). Moreover, 
the spacial sensitivity characteristics of the cable sensor changes along the sensor in 
dependents of the dielectric relaxation behavior of the surrounding material. For this 
reason, the precise determination of soil moisture profiles requires an improved TDR- 
waveform analysis strategy under consideration of the change of the sensitive area along 
the cable sensor. 

A disadvantage of the used FEM with HFSS is the lack of consideration an 
exchange of material between air or water- filled gap and surrounding medium. In 
addition, to achieve sufficient accuracy the numerical simulations with HFSS requires 
a high degree of computational cost in terms of computational time and memory. 
Further, theoretical, numerical, and experimental investigations in conjunction with 
reconstruction algorithms have to point out to what extent the accuracy of water content 
and porosity profiles can be determined in strong dispersive, high lossy materials. 
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